Dimensional Crossover in Spin Diffusion: A manifestation of the Quantum Zeno Effect 
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The Quantum Zeno Effect (QZE) implies that a too frequent (u^ — > oo) observation of a quantum 
system would trap it in its initial state, even though it would be able to evolve to some other state 
if not observed. I n our scheme, interacting spins in a 3-d cubic lattice, "observe" each other with a 
frequency oc -J ' J| + J| -+- Jf /%, where J's are the coupling constants. This leads to a "diffusive" 
spread of a local excitation characterized by the constants oc J^/lo^. Thus, a strongly asymmetric 
interaction (e.g. J y /J x ( z ) 3* 1), would hinder diffusion in the perpendicular directions (D x (z) ~ * 0) 
manifesting the QZE. We show that this effect is present in numerical solutions of simple 2-d systems. 
This reduction in the diffusion kinetics was experimentally observed in paramagnetic compounds 
where the asymmetry of the interaction network manifests through an exchange narrowed linewidth. 
New experimental designs are proposed. 



Quantum dynamics of magnetic excitations in a sys- 
tem of interacting spins at high temperature is an ac- 
tive field of researchfl. From a macroscopic scope, i.e. 
for long times and large wave-lengths, one expects that 
excitations should evolve irreversibly. The resulting hy- 
drodynamic equationso describe the "spin diffusion" . In 
this regime a local excitation decay as (Dt)~ d ^ 2 where 
d is the dimension of the space. From the microscopic 
point of view, however, dynamics is governed by re- 
versible quantum mechanics as long as quantum coher- 
ence is maintained. Although this dynamics seems to 
be "diffusive" already at intermediate times, a careful 
studyu of low dimensional systems can recogniz£.x[uan- 
tum interferences. Its experimental observatiorod con- 
stitutes a fingerprint of the bounded regions where the 
dynamics occurs. Besides™ an experimental realization 
of a "Loschmidt daemon"B, which allows the evolution 
backwards in time, can be achieved by inverting the sign 
of the effective Haniiltonian. In this case, "irreversible" 
interactions are simply those we do not control. Hence, 
reversibility is not totalQ. This effect is stronger than 
what can be inferred from the magnitude of the non- 
inverted terms. In fact, the apparently "diffusive" dy- 
namics of the many body interaction seems to transform 
small residual interactions intcuefficient mechanisms to 
stabilize an irreversible diffusionO. Irreversible effect is of- 
ten amplified even by those interactions we can control. 
This is analogous to the Drude approximation for the 
electrical resistance of an impure metal at low tempera- 
tures. There, the reversible elastic scattering with impu- 
rities with rate (1/T imp ) facilitates the interaction with 
the thermal bath. This bath acts through the uncon- 
trollable electron-phonon interactions or other dephasing 
collisions with rate l/r^. The remarkable consequence is 
that in first approximation the diffusion constant does 
not depend on but it is D oc t> 2 Ti mp , where v is a typi- 
cal velocity for ballistic propagation of the excitation. Its 
lesson is that a reversible interaction can provide an evo- 



lution close enough to diffusion. Then, when actual irre- 
versible processes occur, they stabilizccl the evolution into 
an irreversible diffusion. The observed diffusion constant 
remains the unchanged. This gives the ultimate justifica- 
tion for the stosszahlansatz or assumption of randomness 
after successive collisions in which no memory of the pre- 
vious quantum state is retained. This is equivalent to 
consider a collision with an impurity as a classical mea- 
surement or wave function collapse. This, of course, is an 
approximation which breaks down close to the localized 
regime where interferences play a fundamental role. 

The goal of the many-body techniques is to provide the 
match between the quantum and hydrodynamic regimes 
providing the machinery for the calculation of the diffu- 
sion constant and other macroscopic observables. How- 
ever, because of the implicit approximations, this is some- 
times done at the cost of the physical insights into the 
nature of irreversibility. To understand the dynamics 
of a local excitation, we notice that one spin up, in a 
lattice with all spins down, propagates with a typical 
"ballistic" velocity proportional to the exchange constant 
J. At high temperature, there are as many up spins 
as down ones. Hence, this propagation in the sublat- 
tice of spins down is interrupted by the modification of 
this lattice. This produces a "collision" rate (l/r m b) also 
proportional to J. If we neglect quantum interferences 
{stosszahlansatz) produced by multiple collisions, a dif- 
fusion constant D cx i> 2 T m b oc J is obtained. Again the 
point is that a diffusive behavior is a good approximation 
for the dynamics. It will be stabilized by later irreversible 
interactions. Keeping this in mind, we may consider each 
"collision" as a measurement process. 

Quantum dynamics is strongly modified by recurrent 
measurements. This phenomenon, known as the Quan- 
tum Zeno Effect (QZE)E3, has been applied with partic- 
ular success in quantum opticsEj. In simple words, it 
affirms that if the evolution of a quantum system is ob- 
served too frequently, there will be no evolution to be seen! 
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While this might sound paradoxical from a classical point 
of view, within the quantum logic this is not a paradox at 
all, since a measurement involves a collapse of the wave 
function and the start of a new quantum evolution. For 
short times, the probability of staying in the initial state 
is Pi,i(t) = 1 — (Jt/2U) 2 + . . . , with J being an average 
exchange energy. The lack of a linear term on this expan- 
sion has very important consequences. If the evolution 
in a time t s is interrupted by iV observations, the final 
probability of stay is P M (i = t s ) = [Pis(t s /N)] N -> 1 
when N — > oo. Here, t s /N = defines the decoher- 
ence time. Hence one sees why a fpeouent collapse ( i.e. 
cj0 = 1/t0 » 1) hinders evolutior£3. The result of the 
successively interrupted evolution is shown schematically 
in Fig. 1. 

This article addresses a paradoxical aspect of the di- 
mensional dependence in the dynamics of magnetic exci- 
tations at the light of the QZE. In order to make it more 
obvious, let us firstly discuss the particular limit aa cu- 
bic lattice where the interaction along direction y grows 
(J y /J x (z) ~ * oo)- The diffusion coefficient in that direc- 
tion will also grow: D y oc J y . However, the spreading 
rate in the other directions will decrease D x ,(z) ^ V^y 
The essence of the argument that we are going to de- 
velop is that spins spreading in each xz-plane register the 
evolution in the other parallel planes through the many 
body coupling J y (i.e. they "observe" each other with a 
frequency lo^ — > J y /Ti ), producing the reduction of the 
dynamics within the planes. This is a general behaviour 
also valid for the non-pertuvative limit of J y 
our knowledge this has not been noticed previously. 

The calculation of the dynamics of a local spin polar- 
ization in the high temperature regime can be mapped 
to a system of Fermi particles on a lattice. Up spins are 
identified with particles and down spins with holes. A 
particle is localized at the excited site with the rest of 
the sites being occupied with probability one half. The 
initial state is an incoherent superposition of all the pos- 
sible initial states a few of which are shown in the inset 
of Fig. 1. The particle initially at Oth site starts to 
evolve having a finite probability amplitude to jump into 
empty neighboring sites. Meanwhile, other particles also 
can move in and out of neighboring sites creating a fluc- 
tuating effective potential. If this is approximated by 
a stochastic potential, one obtains an irreversible equa- 
tion satisfying the hydrodynamic limits. The essential 
point is that the correlation times of the potential are the 
same that characterize the dynamics of a particle in its 
fluctuating environment. This leads to a self-consistent 
equation which is the core of classic many-body calcu- 
lations o£-.spin "diffusion" such as that of Blume and 
HubbardEJ for symmetric lattices. However, those cal- 
culations do not show the mind-teasing behavior of the 
asymmetric lattices which are the purpose of this work. 
Let us do a simplified calculation valid for both symmet- 
ric and asymmetric systems, which will clarify the effect 
of dimensional crossover in the dynamics. Consider a d- 
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dimensional (hyper)-cubic lattice of spins 1/2 interacting 
through a Hamiltonian: 
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where S are the usual spin operators with subscripts in- 
dicating spin sites and Jjk — Jp, are nearest neighbors 
interaction parameters depending only upon the direc- 



tion ft along sites j and k at distance a p 



1. For 



a = — i, it describes the Heisenberg model (isotropic ex- 
change), a — defines the XY model, while a — 1 is a 
truncated dipolar Hamiltonian. 

For each initial state, | i),with 0-th site polarized (i.e. 
one of the states shown in Fig. 1), the probability of 
finding the same site polarized in the state (/ | after a 
time i, is 



(f | exp[--Hnt] | i) 



(2) 



A total ensemble averaged probability that a spin initially 
up at position 0-th is still up at time t can be calculated 
summing over all the TV; and Nf possible initial and final 
states: 

N f Ni 



/ i 



(3) 



(Zfj) = 1 is the average number of neighbors along di- 
rection [i with spin down. In the second order term only 
the flip-flop terms produce the exchange of the originally 
polarized 0-th spin with its neighbors. This corresponds 
to the one body dynamics of one up spin in a lattice in 
which the other up spins remain frozen. Higher order 
terms contain the dynamics of those other spins and the 
many body interactions. 

A normalized magnetization can be calculated from 
the spin autocorrelation function as 



M(t) = (S$(t)SZ)/(S$S$) 



2((P(t)) - \). (4) 



This magnitude is experimentally accessible. For short 
times, the mean square displacement of the magneti- 
zation in terms of the nearest neighbor spins correla- 
tion functions: (rfy oc a 2 ^oio^W^oy • The truncated 
quantum dynamics given by Eq. (H) gives : 
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2 t2j.2 /t-2 



(5) 



However, the perturbative expansion of Eq.(||) is not a 
practical way to obtain the long time dynamics but for 
few simple one dimensional systemsEJ. Therefore, we need 



2 



an entirely different approach. First, we notice that for 
infinite lattices the long time evolution is very complex 
and it presents the apparently diffusive behavior (even- 
tually stabilized by interactions with the thermal bath) 
that we want to evaluate. While evolution with ([!]) is not 
an irreversible process, it is close enough to a diffusive 
evolution so that a self-consistent condition is already 
achieved for intermediate times. 

When (P(t)) has decayed substantially, let us say when 
the second term in Eq.(^) is half the first, the environ- 
ment has changed completely and no phase coherence 
with the initial one-body state is retained. This defines 
the dephasing time 



t = mi 



\ 



(6) 



ft is important to note that this functional dependence 
on J ^ does not depend on the amount of decay chosen to 
determine t$. The change in the environment, following 
itself a quantum dynamics according to Eq. (||) , is slower 
than linear at early times, becoming important only at 
about T0. Then, it can be described by a discrete time 
Markovian process leading to a classical random walk. 
At the dephasing time the mean square displacement of 
the magnetization is 



(r 2 ) 



(7) 



At t = t,],, both quantum and classical diffusive regimes 
must coincide. This statement was rigorously pr-e.ved 
in Ref.El, where we used the Keldysh's formalismli-3 to 
achieve a non-perturbative description of the cross-over 
from the quantum to the diffusive regime for a particle 
interacting with a dephasing field. In our model this 
field acts at a typical time with an interaction probabil- 
ity p(t)dt = S(t — T^dt. Then, the coherence with the 
initial state has a survival probability 9{t^ — t). This step 
function is more appropriate to describe the quantum dy- 
namics of the dephasing field than the usual exp[— t/r^]. 
According to Ref.i, within this approximation the self- 
consistent propagation of density excitations (satisfying 
the integral Keldysh equation) requires that both quan- 
tum and Markovian descriptions^ give the same proba- 
bility distribution at t = t^. This condition is equivalent 
to equate (||) and (]?]), from which we obtain the diffusion 
constant for each direction: 



D - x 7 2 / 

2h -V 
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(8) 



This important result contains the paradoxical aspects 
of spin dynamics we discussed in the introductory para- 
graphs. While our procedure has been mainly quali- 
tative, we believe it catches the fundamental phenom- 
ena, and therefore the correct functional dependence. 
For the symmetric three dimensional lattice this gives 



D = J a 2 /(7i2v6), in fair agreement with-tie values cal- 
culated by a number pf. previous authoraLj'Ej and consis- 
tent with simulationst2l in a system of classical spins. 

While numerical solution of the quantum dynamics of 
systems with a large number of spins is a formidable task, 
we can attempt to see signatures of the discussed phe- 
nomena for small two dimensional systems. We consider 
a nine spins system with periodic boundary conditions 
and evaluate the eigenstates of all the spin configura- 
tions to study the dynamics according to Eqs. (|f|) and 
(||) with a = —1/2. In Fig. 2 we show the numerical 
evaluation of (r 2 ) and (r 2 ) as a function of time t for 
different values of J y while keeping J x = 1. . While the 
diffusive regime {{r 2 ) ~ 2D y t) can not be reached in a 
small system, the plot shows that an increase in J y is 
correlated with an increase in the spreading dynamics. 
The paradox that manifest the QZE is that the growth 
of (r 2 }) , showing the dynamics in the perpendicular di- 
rection, is slowed down by the increase of J y . While 
eventually the hydrodynamic limit could be described by 
a diffusion equation in which variables can be separated, 
in the many body Schrodinger equation variables appear 
intimately entangled, leading to the interdependence of 
the diffusion constants of Eq. (^) . 

We want to show that even when an asymmetry of the 
lattice could lead to a faster quantum decay of the po- 
larization for short times, it produces a reduction of the 
effective dimensionality of the lattice where diffusion oc- 
curs which slows down the spreading. For this purpose 
we study the time decay of the local polarization in a 
square lattice. The inset in Fig. 3 shows the magneti- 
zation M(t) in the symmetric lattice up to intermediate 
times (for long times weak mesoscopic beatsa resembling 
to those of 1-d rings would appear). The symmetry of 
the interaction network can be broken by slightly increas- 
ing the coupling along direction y, J y = J + 8 J, while in 
the other direction it is decreased by the same amount 
J x = J — 5 J. The thick line in the main plot shows 
the difference between the local magnetizations calcu- 
lated for the asymmetric and the symmetric networks 
SC = M A {i) - M s (t), with 6 J/ J = 0.1. For very short 
times the spreading in the asymmetric lattice is faster 
than in the symmetric one and SC follows the parabolic 
approximation SC w —{SJ) 2 t 2 /Ti 2 shown by the thin line. 
QZE manifests for intermediate times, slowing down the 
diffusion in the asymmetric case. This compensates the 
fast decay of short times at around t — 1.8%/ J, when 
SC = 0. This crossover from faster decay of M A as com- 
pared with M s to a slower one is a very remarkable result 
of our simple theory consistent with the numerical solu- 
tions. 

A simple experimental test of the spin dynamics is the 
linewidth of a magnetic resonance spectrum. Different 
local environments in a system of non-interacting spins 
produce an inhomogeneously broadened absorption line, 
where each frequency corresponds to spins seeing a dif- 
ferent local field. However, through the flip-flop mecha- 
nism, each spin excitation explores different lattice sites 
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producing-an averaging of the field. Thus, the narrowed 
linewidthEa is: Av oc J °° (S^(t)S^) dt (eventually the in- 
tegral extends up to a cut-off time). Then, the slower the 
dynamics the wider the line. 

It is not usually possible the control of the magnitude 
of the exchange interaction for electronic spins. How- 
ever, this is just the situation observedtll for a family of 
Cu(aa)2 single crystals where aa stands for amino acid. If 
the aa is present in a mixture 50% dextrogyre and 50% 
levogyre (D,L), the Cu 2+ paramagnetic centers occupy 
sites with inversion symmetry in a two-dimensional cou- 
pling network. The Cu atoms are connected by OCO 
bridges, with the different aa residues playing the role of 
separators between layers. In Cu(L-aa)2 crystals (100% 
L-aa), however, a breakdown of this symmetry occurs. 
The amount of asymmetrization of the lattice depends 
on the aa. The very important increase of the electron 
paramagnetic resonance linewidth that is observed when 
going from D,L to L crystals can be interpreted as a sig- 
nature of the discussed QZE. 

Systems with interacting nuclear spins seem quite 
promising for the study of dynamics in asymmetric lat- 
tices. Macro-molecules could be engineered to present a 
sequence of through bond isotropic couplings when stud- 
ied in solution. In crystals, where the dipolar interaction 
is the dominant one, its dependence on the angle 9 be- 
tween the vector connecting dipoles and dipole orienta- 
tion (fixed by an external field) can be used to change 
its magnitude and sign. For example, in a cubic lat- 
tice of spins, by varying the magnetic field in the plane 
[001] from the [110] direction toward one at a magic an- 
gle (9 m = arccos [1/V3]) with the [010] axis, one could 
see how the dynamics changes from a three-dimensional 
behavior (J x = —d/2, J y = —d/2, J z = d), with a sym- 
metric diffusion in the xy-plane, toward that of a two 
dimensional system (J x = —d, J y — 0, J z — d). The 
change from the initial orientation would allow a study 
of a crossover in SC similar to that discussed above. The 
application of these concepts would be even more direct 
for magnetically two-dimensional systems ( J z « 0) where 
J x and J v can be controlled independently. In both cases 
multiple quantum coherences experiments and spin dif- 
fusion pulse sequences] could be used to obtain comple- 
mentary informationta about the asymmetric spreading 
of magnetic polarization. A more conventional experi- 
ment is to apply a magnetic field gradient along one crys- 
tal direction (say y) which, by detuning the resonance fre- 
quencies of nuclei at different xz planes, is equivalent to 
a decrease of the inter-plane coupling. Results obtained 
using this techniqueEj are also consistent with our pre- 
dictions. In summary, we have put the problem of spin 
diffusion under a new perspective that could stimulate 
new series of experiments and calculations to understand 
spin dynamics. 
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Figure 1. Probability of stay for a particle whose 
quantum evolution is interrupted by a measurement pro- 
cess with a period and = 2t^. The coherent evolu- 
tion is approximated up to the quadratic term and the 
return probability is neglected. The increase of the prob- 
ability of stay is a manifestation of the QZE. The evolu- 
tion shown is equivalent to the average evolution of the 
ensemble in the inset. Only 5 of the 16 configurations 
corresponding to a local excitation in the high tempera- 
ture limit are sketched. 

Figure 2. Mean square displacement of magnetization 
as function of time for increasing asymmetrization of the 
2-d lattice with Heisenberg interaction. The full line is 
J y =J x =J, the other lines are the sequence J y /J—2, 3, 4, 
5 and 6. J x is kept constant. Upper curves are displace- 
ments along y, while lower ones are along x. The QZE 
manifests in the slow down of the spreading along x. 

Figure 3. The inset shows the spin autocorrelation 
function in a symmetric 2-d lattice with Heisenberg in- 
teraction as function of time. Thick curve in the main 
frame shows the difference, SC, between the autocorre- 
lation function of the symmetric (J y =J x =J) and asym- 
metric (J v = J + SJ and J x = J — SJ, with SJ/J = 0.1) 
lattices. The thin line is a parabolic approximation. 
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